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Claims 

What is claimed is: 

1 . A method of designing an interbody fusion cage comprising: 
defining operational parameters for the cage; 

defining a macroscopic structural layout for the cage satisfying the 
operational parameters; 

dividing the macroscopic structural layout of the cage into a plurality of 
discreet sub-segments; 

defining a density distribution of the macroscopic structural layout by 
determining a density level for each sub-segment; and 

defining a microscopic structural layout for the cage by assigning pre- 
selected microstructures to the sub-segments in accordance with the density 
level of each sub-segment. 

2. The method of claim 1 wherein the operational parameters 
further comprise stability, porosity, and compliance. 

3. The method of claim 1 wherein the step of defining the 
macroscopic structural layout further comprises executing a topology 
optimization algorithm for the cage based on the operational parameters. 

4. The method of claim 1 wherein the step of defining the 
microscopic structural layout further comprises executing a topology 
optimization algorithm for the cage based on the density distribution. 

5. The method of claim 1 further comprising integrating the 
microscopic structural layout and the macroscopic structural layout to provide 
a designed cage. 

6. The method of claim 1 further comprising manufacturing the 
designed cage with solid free-form fabrication techniques. 

7. The method of claim 1 further comprises categorizing the sub- 
segments into different ranks based on the density level of each sub-segment, 
each rank being defined by a different length scale; and 

homogenizing the microstructure of a particular rank to an upper rank. 

8. A method of designing an interbody fusion cage comprising: 
modeling a mechanical environment in which the cage is to be used; 
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performing a macroscopic layout topology optimization process to 
define a global layout topology solution for the cage; 

performing a density distribution process on the global layout topology 
solution to define a density distribution of the global layout topology; 
5 performing a microscopic layout topology optimization process to 

define a microstructure topology solution according to the density distribution; 
and 

integrating the macroscopic and microscopic topology solutions. 

9. The method of claim 8 wherein the step of modeling the 
10 mechanical environment further comprises creating a finite element model to 

simulate the mechanical environment. 

10. The method of claim 8 wherein the step of performing a 
microscopic layout topology optimization process to define a microstructure 
topology solution according to the density distribution further comprises 

15 representing microstructures of sub-segments of the cage with intermediate 
density values from the global topology optimization solution, the 
microstructures being defined at different length scales than the scale of the 
global optimization solution. 

11. The method of claim 8 wherein the step of performing a 
20 microscopic layout topology optimization process to define a microstructure 

topology solution according to the density distribution further comprises: 

employing a low resolution mesh with a homogeneous density 
as a first initial guess in the optimization process; 

solving an optimization problem for the low resolution mesh to 
25 yield a rough solution; 

meshing the rough solution at a finer resolution to yield a finer 
resolution solution; 

using the finer resolution solution as a subsequent initial guess 
in the optimization problem; and 
30 repeating the above steps until a finest resolution solution is 

provided. 
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12. The method of claim 8 wherein the step of performing a 
microscopic layout topology optimization process to define a microstructure 
topology solution according to the density distribution further comprises: 

defining one of stiffness and porosity as an objective function of the 
5 cage and the other of stiffness and porosity as a design variable; 

inputting material constraints for the cage based on the objective 
function; 

inputting a low resolution mesh of the microstructure topology based on 
the objective function as an initial guess design; 
10 homogenizing the material constraints and calculating sensitivity 

functions for the initial guess design; 

determining if a convergence of the objective function is achieved by 
the initial guess design; 

if convergence is not achieved: 
15 solving a topology optimization problem in the optimization 

process using a method of moving asymptotes; 

updating the design variable in the initial guess design; and 
repeating until convergence of the objective function is 

achieved; 

20 once convergence is achieved, filtering the initial guess design; 

determining whether a finest mesh size has been achieved within the 
initial guess design; 

if the finest mesh size is not achieved, updating the material constraints 
and refining the low resolution mesh to yield an updated initial guess; 
25 applying image processing filtering techniques to the updated initial 

guess; 

providing the image processing filtered initial guess for a smaller mesh 
size as a new initial guess design; and 

repeating the above steps until the finest mesh size is achieved. 
30 13. The method of claim 12 further comprising performing heuristic 

random rounding filtering on the finest mesh size design to yield a final design 
layout. 
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14. The method of claim 12 wherein the image processing 
technique eliminates checkerboard density patterns. 

15. The method of claim 14 wherein the image processing 
technique employs a Gaussian smoothing filter and a connectivity filter to 
smooth sub-segment densities within the cage relative to surrounding sub- 
segments. 

16. The method of claim 12 wherein the stiffness objective function 
further comprises: 

Objective junction: 



Min 

E**™* ,d h d 2 ,d 3 




^.bone eff - tissue eff 



bone eff 



cl 



£ bone eff £, scaffold eff ' 



bone eff 



where n = 1-9. 
Constraints: 
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d,,d 2 .d 3 < 900 
di,d 2 ,d 3 > 300 um 
> % porosity, 

Vtotal 

^scaffold ^ Emin, 

gscaffold ^ jE;maX) 

wherein the material constraints include: 

£ sca/foW as the scaffold base material Young's modulus; 
di, d2, and d 3 as three cylindrical diameters of pores within the 
microstructure; 

d™* 0ff as the effective stiffness of the target bone; 
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tissue eff 



as the regenerate tissue effective stiffness; and 



cscaffoideff^ the scaffoId effective stiffness. 
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1 7. The method of claim 1 2 wherein the porosity objective function 
further comprises: 

Objective Junction: 

Max Vpore , 

,d h d 2 ,d 3 viotal 

Constraints: 

x \ C .^° ne e ^ < ^^ ssue e ^ < X 2 C .^ one e ^ 

where i = 1-9; x 2 > jcj , 
^ ^bone eff < ^scaffold eff < ^ C i bone eff 

where i = 1-9; ^>^, 
dj t d 2 ,d 3 < 900 M-m, 

dj t d 2 ,d 3 > 300 nm, 
^scaffold > £;min, 
^scaffold - £max, 

wherein cti, a 2 , pi, and p 2 are scaling factors used to bound the cage 
5 and regenerate tissue effective stiffness and the material constraints include: 

£ scaffoW as the scaffold base material Young's modulus; 

du d2, and d 3 as three cylindrical diameters of pores within the 
microstructure; 

Qtoneen as the effective stiffness of the target bone; 
1 0 gissue eff as the re g enera te tissue effective stiffness; and 

(jscaffoideff^ thQ 3^01^ effective stiffness. 

18. A method of designing an interbody fusion cage comprising: 
generating a global density distribution under physiologic loading for 
the cage using a global topology optimization algorithm including: 
1 5 a stability constraint that limits total displacement of the cage at 

a desired surface to be less than a target value; and 

a total porosity constraint that ensures desired biofactor 
delivering ability and compliance; 

segmenting the global density distribution architecture into a plurality of 
20 regions, each region having a material phase selected from: 
a low porosity solid phase; 
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a high porosity solid phase; and 
a voided phase; and 
defining a porous microstructure for the cage by generating periodic 
microstructures for the regions having the high porosity solid phase and low 
5 porosity solid phase using a microstructure topology optimization method. 

1 9. The method of claim 1 8 wherein: 

the low porosity solid phase regions are about 45%~55% solid; 
the high porosity solid phase regions are about 20%~35% solid; 

and 

10 the voided phase regions are about 0% solid. 

20. The method of claim 18 wherein the step of defining a porous 
microstructure for the cage further comprises achieving Hashin-Shtrikman 
stiffness bounds for porous isotropic materials. 
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